Astron. Nachr./AN xxx (xxxx) x, xxx-xxx 



Simulating large-scale structure formation with magnetic fields 



Klaus Dolag 

Max-Planck-Institut fur Astrophysik, P.O. Box 1317, D-85741 Garching, Germany 



Received; accepted; published online 



Abstract. In the past, different works based on numerical simulations have been presented to explain magnetic fields (MFs) 
in the large scale structure and within galaxy clusters. In this review, I will summarize the main findings obtained by different 
authors and - even if many details are still unclear - 1 will try to construct a consistent picture of our interpretation of large- 
scale magnetic fields based on numerical effort. I will also sketch how this is related to our understanding of radio emission 
and summarize some arguments where our theoretical understanding has to be improved to match the observations. 
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1. Introduction 

There are several different approaches to study the build- 
up of MFs in the intergalactic medium. Many papers 
consider MHD sim ulations of cloud-wind interaction (see 
Oreg on etalJ 200(1 and references there in) or are si mulat- 
ing the rise of relic radio bubbles (see Uones & De Younj 
2005, Reynol d! et all20 05 and references herein). Such work 
is usually focused more on the relevance of the local MFs 
within these processes, whereas in this review I will con- 
centrate more on the discussion of MHD simulations which 
aim at understanding the build-up of the cluster MF and its 
possible origin. So far, firm evidence for the presence of ex- 
tended MFs has been found on l y in ga laxy c lusters. For recent 
review s see ICarilli & TavloJ (|2002) and iGovoni & Fe retti 
(2004). Recently there are observational claims of a detec- 
tion of substantial MFs also within the large scale structure 
(LSS), see contribution by Kronberg. The origin of the MF 
in galaxy clusters is still under debate. The variety of pos- 
sible contributors ranges from primordial fields, battery and 
dynamo fields over all classes of astrophysical objects which 
contribute with their ejecta. The later possibility is supported 
by the observation of the metal enrichment of the intracluster 
medium (ICM), which has to be originated as ejecta of astro- 
physical objects. In addition, the MFs produced by all these 
contributors will be compressed and amplified by the process 
of structure formation. The exact amount of this amplification 
and the resulting MF filling factor will depend on place and 
time at which the contributer is thought to be most efficient. 
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2. Possible origins of MFs within the LSS 

In the following I will describe the three main classes of mod- 
els for the origin of cosmological MF within the LSS. 

In the first one, MFs are assumed to be produced 'lo- 
cally' at relatively low redshifts (z ~ 2 — 3) by galactic (e.g. 
Volk&Atovan|2000|) or AGN ejecta (e.g. Furlanet to & Loebl 
2001). One of the main arguments in favor of this model is 
that the high metallicity observed in the ICM suggests a sig- 
nificant enrichment driven by galactic winds or AGNs in the 
past. Winds and jets should carry MFs together with the pro- 
cessed matter. While it has been shown that winds from or- 
dinary galaxies give rise to MFs which are far weaker than 
those observed in galaxy clusters, MFs produced by the ejecta 
of starburst galaxies can be as large as 0.1 /iG. Clearly, this 
class of models predicts that MFs are mainly concentrated in 
galaxy clusters. Note that, if the magnetic pollution happens 
early enough (around z ~ 3), these fields will not only be 
amplified by the adiabatic compression of the proto-cluster 
region, but also by shear flows, turbulent motions, and merg- 
ing events during the formation of galaxy clusters. 

In the second class of models, the MF seeds are assumed 
to be produced at higher redshifts, before galaxy clusters 
form gravitationally bound systems. Although the strength of 
the seed field is expected to be considerably smaller than in 
the previous scenario, the adiabatic compression of the gas 
and the shear flows driven by the accretion of structures can 
give rise to a considerable amplification of the MFs. Several 
mechanisms have been proposed to explain the origin of MF 
seeds at high redshift. Some of them are similar to those dis- 
cussed above, differing only in the time at which the mag- 
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Fig. 1. Example of MF (y-component) within an evolving KH insta- 
bility (inlay) and its time evolution. Taken from Birk et al. ( 1999). 



netic pollution is assumed to take place. In the present class 
of models the MF seeds are supposed to be expelled by an 
early population of dwarf starburst galaxies or by AGN a t 
a redshift between 4 and 6 feronberg. Lesch & Hopplfl 999). 
allowing them to magnetize a large fraction of the volume. 
Alternative models invoke processes that took place in the 
early universe. Indeed, the ubiquity of MFs in the universe 
suggests that they may have a cosmological origin. In general, 
all 'high-z models' predict MF seeds filling the entire vol- 
ume of the universe. However, the assumed coherence length 
of the field crucially depends on the details of the models. 
While scenarios based on phase transitions give rise to coher- 
ence lengths which are so small that the corresponding fields 
have probably been dissipated, MFs generated at neutrino or 
photon decoupling have much higher chances to survive un- 
til the present time. Another (speculative) possibility is that 
the seed field was produced during inflation. In this case, the 
coherence length can be as large as the Hubble radius. See 
Grasso & Rubinsteii] fcOOlh for a recent review. 

The third scenario assumes that the MF se eds were pro- 
duced by the so-called Biermann bat tery effect ( Kulsrud et al.l 
ll997lRvu. Kan? & Biermannlll998h . The idea here is that 
merger shocks related to the hierarchical structure formation 
process give rise to small thermionic electric currents which, 
in turn, may generate MFs. The battery process has the attrac- 
tive feature of being independent of unknown physics at high 
redshift. Its drawback is that, due to the large conductivity of 
the IGM, it can give rise to at most very tiny MFs, of order 
10 -21 G. One therefore needs to invoke a subsequent turbu- 
lent dynamo to boost the field strength to the observed level. 
Such a turbulent amplification, however, cannot be simulated 
numerically yet, making it quite difficult to predict how it 
would proceed in a realistic environment. It is clear that one 
expects the level of turbulence to be strongly dependent on 
the environment, and that it should mostly appear in high- 
density regions like collapsed objects. While energetic events 
such as mergers of galaxy clusters can be easily considered 
to drive the required levels of turbulence, this is harder to 
understand in relatively quiet regions like filaments. Lacking 
a theoretical understanding of the turbulent amplification, it 
is therefore not straightforward to relate the very weak seed 








Fig. 2. The gas velocity field in a slice through the central Mpc of 
a cluster simulation after subtracting the global mean bulk veloc- 
ity of the cluster. The left panel shows a run with the original SPH 
viscosity, the right panel for a low-viscosity scheme. The underly- 
ing color maps represent the ratio of turbulent kinetic energy and 
total kinetic energy content of particles, inferred after substracting 
the local mean velocity, as described in lDolag et alJ ( 2005) 



fields produced by the battery process to the MFs observed 
today. Attempts to construct such models based on combin- 
ing numerical and analytical computations have not been re- 
ported to successful reproduce the observed scaling relations 
of MFs in galaxy clusters so far. 

3. Local amplification of MFs 

A very basic process of amplification for MFs is related to the 
Kelvin-Helmholtz (KH) instabilities driven by shear flows, 
which are commo n within the dynamic s of the structure for- 
mation process. iBirk. Wiechen & Qttol £l999) performed a 
detailed study of such an amplification within the environ- 
ment of galactic outflows (see Fig. [T} in starburst galaxies 
where the KH timescale should be ~ 4 x 10 5 years. Using 
a Cartesian resistive MHD code they found that the obtained 
amplification factor for the MF mainly depends on the initial 
ratio of magnetic to kinetic energy and only mildly depends 
on the assumed resistivity. They reached the conclusion that 
such a process could indeed explain the significantly higher 
MF observed in the starburst galaxy halo compared to what is 
expected from the MFs observed within galactic disks. When 
applied to a cluster core environment, the KH timescale turns 
out to be 10 7 years, making it an interesting process for fur- 
ther amplifying weak MFs. 

In recent high-resolution smoothed particle hydrody- 
namic (SPH) simulations of galaxy clusters within a cos- 
mological environment, using a novel sche me to treat ar - 
tificial viscosity w ithin the Gadget! code (Sprin gell2QQ5l) . 
iDolag et all (120051) demonstrated how such shear flows, 
which are quite common in the process of cosmic structure 
formation, drive fluid instabilities and increase the turbulence 
level within the ICM to a significant level (see Fig.|3. It was 
also shown that the artificial viscosity used in standard SPH 
simulations might significantly suppress such fluid instabili- 
ties. 

In earlier work, extensive MHD simulations of single 
merging events performed using the Eulerian code ZEUS 
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Fig. 3. The evolution of (logarithm of) gas density (left column), 
gas temperature (central column ), and (logarithm of) magnetic pres- 
sure (right column) in two-dimensional slices taken through the 
cluster core in the plane of the merger. Each row refers to different 
epochs: t = (i.e. the time of core coincidence), t = 1.3, t = 3.4, 
and t = 5.0 Gyrs, from top to dow n. Each panel is 3.75 x 3.75 Mpc. 
Taken from Roettige retafl JlQQQh . 



JStone & Normanlh 992albl) demonstrated the MF amplifica - 
tion in such merger events (Roettiger, Stone & Burnsfll999). 
In particular they found that first, the field becomes quite fil- 
amentary as a result of stretching and compression caused by 
shocks and bulk flows during infall, but only a minimal am- 
plification occurs. Second, field amplification is more rapid, 
particularly in localized regions, as the bulk flow is replaced 
by turbulent motions (e.g., eddies), see Fig [3] The total MF 
energy is found to increase by nearly a factor of three with 
respect to a non-merging cluster. In localized regions (associ- 
ated with high vorticity), the magnetic energy can increase by 
a factor of 20 or more. A power spectrum analysis of the mag- 
netic energy showed that the amplification is largely confined 
to scales comparable to or smaller than the cluster cores: this 
indicates that the core dimensions define the injection scale. 
It is worth to notice that, due to the lack of resolution, the 
previous results can be considered a lower limit on the to- 
tal amplification. Furthermore, it is quite likely that a galaxy 
cluster undergoes more than one such an event during its for- 
mation process, and that also the accretion of smaller haloes 
injects turbulent motions within the ICM: consequently the 
MF amplification within galaxy clusters will be even higher. 




1.10 & 



0.4 0.6 
R (Mpc) 



0.80 
1.0 



Fig. 4. In the left panel we show the simulated A3 667 X-ray sur- 
face brightness and radio data. Contours represent the X-ray surface 
brightness. The grayscale represents the synthetic radio emission at 
1.4 GHz. The image is 3.15 x 3.85 Mpc. Dashed and dotted lines 
refer to the location of the radio spectral index (a\\g) profiles dis- 
played in the right panel. Taken from Roet tiger et alJ dl999l) . 



A detailed discussion of the amplification of MF in a cluster 
environment, usin g various simulations of driven turbulence , 
can be found in |Subramanian. Shuk urov & Haugenl (|2005), 
where it is shown that reasonable strength and length scales 
for galaxy clusters can be obtained by turbulent processes. 



4. The quest for radio relics 

By extending such merger simulations to reproduce 
the observed X-ray properties of A3 667 and adding a 
model for in situ re- accelerati on of relativistic particles, 
Roet tiger. Burns & Stond |l999) were able to reproduce the 
main features of the extended peripheral radio emission (the 
so-called radio relics) observed in A3 667. In their models 
they injected relativistic electrons with a power-law spec- 
trum, where the power-law index 7 = 3/(r — 1) + 1 is related 
to the gas compression rate r at the shock. They also relate 
the age of the radio plasma t a to the distance d = Kv s t a , us- 
ing a weak-field/high-diffusion limit k = 1. Having effective 
shock velocities v s « 700 — lOOOkm/s and aging the syn- 
chrotro n spectrum using the formalism by Mv ers & Spa ngler 
(119851) they have been able to reproduce the observed distri- 
bution of spectral index for a MF of w 0.6/iG at the position 
of the radio relic, see Fig.|4] Since such configurations seem 
to be quite common in galaxy clusters, naturally the ques- 
tion arises, why not all clusters show such a peripheral radio 
emission. One possible explanation is that such shock struc- 
tures are relatively short lived compared to the merger event 
itself. Also, it is necessary to have the presence of a large- 
scale MF. It could further be that only massive clusters can 
provide enough MF and strong enough merger events to trig- 
ger such a peripheral emission. 

To overcome this problem, EnBlin & Briiggen (2002) 
proposed that such radio relics could be made by pre-existing 
fossil radio plasma illuminated via shock waves initiated 
by merger events. In their work they evolved the electron 
spectrum for the tracer particles, representing the fossil ra- 
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Fig. 5. Evolution of the gas density (top) and MF energy density 
(bottom) for a shock interacting with fossil radio plasma. Taken 
from lEnBlin & Briiggenl <|2002 ). 

dio plasma, following EnBlin & Gonal-Krish nal fcOOlh which 
take into account synchrotron, inverse Compton and adiabatic 
energy losses and gains. Their simulation, using the ZEUS 
code, follows the evolution of a sphere of tracer particles hit 
by a shock front (see Fig.|5J). Such a configuration nicely re- 
produces the filamentary radio emission and toroidal struc- 
tures, as observed in many cases. These simulations also pre- 
dict the MFs to be mostly aligned with the direction of the 
filaments, as suggested by observational data. 

This idea was investigated in a more realistic modelling 
bv iHoeft. Brliggen & Yenesl (12004 . using an SPH Gadget 
( Springel, Yoshida & White 2001) simulation of a merging 
galaxy cluster within a cosmological environment. Such a 
simulation showed that the probability for a shock wave to 
flare the radio plasma is highly suppressed in the central re- 
gions of galaxy clusters compared to the peripheral ones, 
where they found illuminated structures up to Mpc in size 
(see Fig.|6l>. The reason for this is that first at the center the 
radio plasma ages much faster due to its higher pressure (and 
the losses coming in together with the higher MF), and sec- 
ond that the compression ratio of the shock wave is much 
higher in the low-density peripheral regions. It is worth to no- 
tice that a necessary condition to form such relics is that the 
initial state of the fossil radio plasma is characterized by a ra- 
tio of Pb / Pgas which is as low as 1 per cent to allow shocks 
to revive « 1-Gyr old radio g hosts. It is also important to 
mention that Hoeft et al1(E()0 4) find high probability of radio 
emission outside of the shocks, related to drained gas flows, 
induced by the merger events, which transport material from 
the outskirts towards the higher density regions. Thereby in 
some cases the adiabatic compression seems to be enough to 
revive the fossil radio plasma. 

5. Shocks in cosmological simulations 

Cosmological shocks, mainly the accretion shocks on cosmo- 
logical objects like galaxy clusters and filaments, are much 
more frequent than the ones produced by individual merger 
events. They also c an produce MFs by the so-called Bier- 
mann battery effect (Kulsrud et al. 1997; iRvu etal.lll998h on 
which a subsequent turbulent may operate (see section |2|>. 




Fig. 6. The projected 'potential' radio luminosities for 1.13 Gyr old 
radio plasma, where Pb/P^&s = 0.01. For comparison the bolo- 
metric surface X-ray luminosity is given. The total bolometric X- 
ray luminosity of the cluster is 2 x 10 44 erg s _1 and the emission- 
weighted temperature is 3 keV. Taken from H oeftetalJ foOOA . 



In such a scenario, the MF is strongly correlated with the 
large-scale structure (see Fig. Q). This means that in such a 
case, the MF within the filamentary structure could be even 
slightly higher than its equipartition value without violating 
the (weak) up per limits of rota tion measure of quasars, as 
pointed out bv lRvu et all Jl998l) . It is worth to point out that 
the arguments for the turbulent dynamo action, which could 
a mplify the battery see d fields up to /iG level, as presented 
in iKulsrud et alJ (Il997h . refer explicitly to regions about to 
collapse into galaxies. It has still to be proven that such argu- 
ments hold within proto-clusters or even cosmological struc- 
tures, like sheets and filaments. In general, the time evolution 
of the MF as predicted by t hese simulations saturates around 
z ~ 3 ( Kulsru d et alJll997|) and leads to a relatively uniform 
MF strength on scales of tens of Mpc within the LSS around 
galaxy clusters (see Fig.0). Note that so far there is no com- 
parison of synthetic rotation measurements obtained by the 
MFs predicted by up-scaling the battery fields, with obser- 
vations on scales of galaxy clusters. This might be partially 
motivated by the lack of resolution in simulated clusters. 

In addition these shocks act as place for acceleration 
of cosmic rays (CR) which then will be accreted into the 
LSS, specially withi n galaxy clusters. Us ing the COSMOR 
(Miniati 2001) code, Mini ati etalJ J200lh followed primary 
ions and electrons (injected and accelerated by diffuse cosmic 
shocks) and secondary electrons and positrons (produced in 
p-p inelastic collisions of CR ions with thermal ICM nuclei) 
within a cosmological simulation. Under the assumption that 
the MF produced by the battery effect reflects a fair represen- 
tation of the true distribution of relative MF strengths within 
the LSS, they were able to predict the central radio emission 
(radio halos) - mainly produced by secondary CRs - as well as 
the peripheral radio emission (radio relics) - mainly produced 
by primary CRs - in a self-consistent treatment (see Fig. 
|9|): the resulting morphology, polarization and spectral index 
match the observed properties. However, one has to notice 
that the extrapolation of these simulations (on group scale) to 



6 COSMOLOGICAL MHD SIMULATIONS 




-20 20 40 

x/Mpc 



Fig. 7. A two-dimens i onal c ut through a cosmological box simu- 
lated by iMiniati et alJ fcOOlh . and following the evolution of bat- 
tery fields. Shown is (the logarithm of) the up- scaled MF strength in 
Gauss taken at z = 0. Taken from lSigl. Miniati & En ss fflirJ J20041 

the observed data (i.e. on massive cluster scale) might not be 
straightforward (see Fig.[OJ>. It is further worth to notice that 
even with a proton injection ratio R P j n of 1Q~ 2 (which is de- 
rived from observations, see Tab.5 in Min iati et all2 001). the 
predicted luminosities for the primary emission is still signif- 
icantly higher than for the secondary one (see Fig.[T3t. which 
seems not to be confirmed by the observations. Although the 
results for the secondary model are i n agreement with previ - 
ous work on massive galaxy clusters (Dolag & EnBlin 2000), 
we notice that generally in the secondary model radio haloes 
are predicted for every cluster, which is in contradiction to 
the observations (see Fig.[T3l. 

6. Cosmological MHD simulations 

Using GrapeMSPH JDolag. Bartelmann & Le"scri ri999) and 
assuming that a small initial magnetic seed field exists 
before structure formation (see section |3, the first self- 
consistent simulations which follow the MF amplifica- 
tion during the formation of galaxy cluster s with cos- 
mological environment have been perf ormed foolag et alJ 
1999 Dolag, Bartelmann & Lesch 2002). These runs were 
able to demonstrate that the contribution to the amplification 
of MFs by shear flows (and by its induced turbulence) is sig- 
nificant (see Fig.fTol). Therefore for the first time a consistent 
picture of the MF in galaxy clusters could be constructed: 
the amplification predicted by the simulations was capable 
to link the predicted strength of the seed MFs (see section 
and references therein) at high redshift (z « 3 and higher) to 
the today observed MF strength in galaxy clusters. Further- 
more the simulations predicted that the final structure of the 
MF in galaxy cluster reflects the process of structure forma- 
tion, and no memory on the initial MF configuration survives: 
this relaxes the constraints on models for seed MFs. In gen- 
eral such models predict a MF profile similar to the density 
profile. Thereby the predicted rotation measure (RM) pro- 
file ag rees with the observed one (see Fig. [Tit . Dolag et al. 
(l200ll) found a quasi linear correlation between two observ- 
ables, namely the X-ray surface brightness and the RM r.m.s.. 
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Fig. 8. Slice through the centre of the simulated box, hosting a 
massive galaxy cluster. The simulation, made using the AMR code 
FLASH, follows the evolution of a weak magnetic seed field. Shown 
is the logarithm of the MF strength in Gauss, as measured at z = 
0.5. Taken from lBruggen et alJ J200A . 



This result is now confirmed to hold over several orders of 
magnitude in a collection of observational data existing in 
the literature (see Fig fT2b. Extendin g Gadget! to follow the 
full set of ideal MHD equations, iDolag et alJ (12004 [2005) 
performed several realizations of a cosmological volume con- 
firming t he previous findings e ven at much higher resolution. 
Recently, Brii ggen et alJ (120051) performed a simulation of the 
formation of a single galaxy cluster in a cosmological frame- 
work, using a passive MHD solver implemented into FLASH. 
Thereby they nicely confirmed all previous results based on 
the SPH codes, using their adaptive mesh refinement (AMR) 
simulation code (see Fig.©. 

Another interesting quantity to look at is the slope a of 
the MF power spectrum (oc k~ a , with k being the wave vec- 
tor). Wit hin galaxy clusters a is predicted by the SP H simu- 
lations iDolag et alJl2002lRordorf. Grasso & Dolagll2004l) to 
be slightly lower, but still very close to 11/3, which is the ex- 
pected value for a K olmogorov-like spectr um (in 3D). The 
AMR simulation by Briig gen et alJ (120051) nearly perfectly 
matches the Kolmogorov slope. Fig. [l4| shows a collection of 
predicted and observed slopes for the MF spectrum. Note that 
this can only be a crude comparison, due to the limitations in 
both observations and simulations. But in general the range 
of slopes predicted by the simulations seems to be consistent 
with the slopes inferred from observations of real clusters. It 
is worth to notice that the numerical results suggest that the 
spread of the slopes might reflect the dynamical stage of the 
system, indicated by the inlay, showing the two mass accre- 
tion histories. Here the smoothly accreting clusters show the 
smallest slope (indicating the lack of large-scale power). On 
the contrary, the clusters having a quite violent history, with 
several major mergers indicated by the small arrows, show a 
very steep spectra. Two of them happened in the last 3 Gyrs, 
thereby inducing an excess of power on large scales. 

Further support for strong MF amplification during the 
process of structure formation comes from the application of 
the Zel'dovich approximation to follow the MHD equations 
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Fig. 9. Synchrotron power at 1.4 GHz from secondary electrons 
(left panel) and primary electrons (right panel). Note that the val- 
ues for the luminosities for primary electrons should be s caled with 
the ele ctron to proton injection ratio R e / P . Taken from Miniati et al. 

teooih . 

durin g the gravitational coll apse feruni. Maartens &Ts agas 
2003; iKing & Colesl |2005). Such works showed super- 
adiabatic amplification due to the anisotropy of the collapse 
of the LSS within the cold dark matter paradigm. 

A novel aspect of including MF pressure into LSS simu- 
lation - even if in a simplified way - is investigated in an on- 
going project (see the contribution by Gazzola), which aims 
at identifying the point for which such non thermal pressure 
support starts to significantly modify the structure formation. 

7. Conclusive remarks 

It seems that within the last years a probably con- 
sistent picture of MFs arises from numerical works 
and observations. Supported by simulations of individual 
events/environments like shear flows, shock/bubble interac- 
tions or turbulence/merging events, a super-adiabatic ampli- 
fication of MF is predicted. It is worth to notice that this 
common finding is obtained by using a variety of different 
codes, based on different numerical schemes. Further sup- 
port for such super-adiabatic amplifications comes from an- 
alytical estimates of anisotropic collapse making use of the 
Zel'dovich approximation. When applied in fully consistent 
cosmological simulations, various observational aspects are 
reproduced (see Figs. [H] and [12|> and the resulting MF am- 
plification reaches a level sufficient to link models predicting 
seed MFs at high redshift with the MFs observed in galaxy 
clusters today. It is important to mention that all simulations 
show that this effect increases when the resolution is im- 
proved, and therefore all the numbers have to be taken as 
lower limits of possible amplification. 

Note that on the other hand there are significant differ- 
ences for the predictions of the MF structures coming from 
different models of seed MFs. In particular there are main dif- 
ferences between the up-scaled, cosm ological battery fields 
(Min iati et al1l200ll Isigl et al.l l2004) and the MF predicted 
from high- resolution simulations of galax y clusters using ei- 
ther AMR feriiggen et alJl2005h or SPH (bolag etalJ ll999. 
120021 12005). In the latter case it is possible to follow the am- 
plification of seed fields within the turbulent ICM in more 
detail. A good visual impression can be obtained by compar- 
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Fig. 10. The MF strength as a function of baryonic overdensity. The 
long-dashed line shows the expectation for a purely adiabatic evo- 
lution, the solid line gives the mean fiel d strength at a give n over- 
density within a cosmological simulation (Dolag et al. 2005). While 
the latter is close to the adiabatic value in underdense regions, there 
is a significant inductive amplification in clusters due to shear flows 
and turbulence, subject however to saturation in the cluster cores. At 
any given density, a large fraction of particles remains close to the 
adiabatic expectation, as shown by the dotted line, which gives the 
median of the distribution at each density. 

ing the regions filled with high MFs shown in Figs.0and[8l 
It is clear that the high MF regions for the battery fields are 
predicted to be much more extended, leading to a flat profile 
around the forming structure, whereas for the turbulent am- 
plified MFs the clusters show a much more peaked MF dis- 
tribution. Part of this difference originates from the physical 
model behind, as the cosmological shocks are much stronger 
outside the clusters than inside. Somewhat more unclear is 
what the contribution of the different numerical resolutions of 
these simulations to these discrepancies is. Calibrating such 
simulations using the MF measured only in the high-density 
regions of galaxy clusters makes it crucial to perform a more 
detailed comparison with all available observations. Note that 
extrapolating the predictions of the simulations into lower 
density regions, where no strong observational constraints ex- 
ist, will even amplify the differences in the predictions for the 
MF structure between the different simulations. 

Moreover, one has to keep in mind that depending on the 
ICM resistivity the MF could be suffering decay, which so far 
is neglected in all the simulations discussed before. Further- 
more a clear lack of the present simulations is that they do 
not include the creation of MF by all the feedback processes 
happening within LSS (like radio bubbles inflated by AGNs, 
galactic winds, etc.): this might alter the MF prediction if 
their contribution turns out to be significant. Also all the sim- 
ulations done so far neglect radiative losses: if included, this 
would lead to a significant increase of the density in the cen- 
tral part of clusters and thereby to a further MF amplifica- 
tion in these regions. Finally, there is an increasing number 
of arguments suggesting that instabilities and turbulence on 
very small scales can amplify the MFs in a relatively short 
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Fig. 11. Comparison of RMs from the simulation with observations 
for Abell clusters, as a function of distance to the closest cluster. 
Smooth lines represent the median values of |RM| produced by sim- 
ulated clusters with masses above 5 x 10 14 M and 3 x 10 14 M . 
The broken line represents the median of combined data taken 
from the independent samples presented in IKim et al] Jl99l and 
Clar ke et all 1200 ll) . We also includ e data (diamonds) for the three 
elongated sources observed in A 1 19 iFeretti et _ajllggd) , and for th e 
elongated source observed in the Coma cluster iFeretti et all 19951) . 



timescale reaching the observed jiG level. However it is still 
unclear how such very small scale fields can be ordered on 
large (up to hundreds of kpc) scales as observed. 

Concerning Radio haloes and relics the picture is only 
partially consistent ( for a more detailed discussion of pri- 
mary and secondary models see Bru nettil2004l and references 
therein). Radio relics seem to be most likely related to strong 
shocks produced by major merging events and therefore pro- 
duced by direct re- acceleration of CRs, the so-called primary 
models. Although some of the observed features like mor- 
phology, polarization and position with respect to the clus- 
ter centre can be resonably well reproduced, there might be 
still some puzzles to solve. On one hand, direct accelera- 
tion of CRs in shocks seems to overestimate the abundance 
and maybe the luminosity of radio relics; on the other hand 
simulations which illuminate fossil radio plasma can pro- 
duce reasonable relics only starting from a small range of 
parameter settings. A similar situation arises for modelling 
the central radio emission of galaxy clusters by secondary 
models. On one hand, the total luminosity seems to be repro- 
duced using reasonable assumptions and also the steep ob- 
served correlation between cluster temperature/mass and ra- 
dio power seems to be reproduced quite well. But such mod- 
els suffer from two drawbacks. The first one is that in the 
framework of such models every massive cluster produces 
a powerful radio halo, but this is not confirmed by observa- 
tions (see Fig. [T3t . The other problems is that the detailed 
radio properties are not reproduced. Firstly, the profile of ra- 
dio emission in most cases is too steep, so that these mod- 
els cannot reproduce the size of the observed radio halos in 
almo st all cases. Secondly, the observed spectral steepening 
(e.g. lGiacintucci et al 12 005) also cannot be reproduced. Note 
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Fig. 12. Correlation of rotation measurements with X-ray surface 
brightness for different galaxy clusters. The data points are obser- 
vations collected from the literature. The colors represent the virial 
cluster temperature taken from litera ture. This is an updated version 
of the correlation presented in Dola g etalJfoOOll) . 



that there is also no indication from observations that clus- 
ters showing radio emission contain higher MFs than the ones 
without observable extended, diffuse radio emission. On the 
contrary, the cluster A2142 has a MF strength similar to the 
Coma cluster (see Fig. [T2l» . but the upper limit on its radio 
emission is at least two orders of magnitude below the value 
expected from the correlation (see Fig. [T3t . Note that both 
clusters are merging systems characterized by the presence 
of two central cD galaxies. This indicates that there should 
be further processes involved or additional conditions to ful- 
fill to produce radio emission. It is worth to notice that recent 
models, based on turbulent acceleration, seem to overcome 
this problem (see[K uo. Hw ang & Irjl2QQ3t iBrunetti & Blasi 
2005; Cassano & Brunetti 2005 and references therein). Note 
that lCassano & Brunettil (l2005h predicted the probability for 
a galaxy cluster to show giant radio halo to be an increasing 
function of cluster mass and reproduced the observed fraction 
of « 30% for massive galaxy clusters. 

There is a big challenge for the next generation of cos- 
mological MHD simulations. Simulations are quite close to 
having the resolution necessary to properly describe the MF 
components down to the observed scales, but on the other 
hand this means to resolve galaxies inside the simulations as 
well. Note that following the dynamics of such structures for 
which the MF might dominate the evolution, is a real chal- 
lenge within LSS simulations. However, we expect that - if 
succeeding in overcoming these limitations - the dynamical 
impact of the MF on regions like the cooling flows at the cen- 
tres of galaxy clusters will be significant and will contribute 
to solve these outstanding puzzles. 
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Fig. 13. Total power of radio halos observed at 1.4 MHz vs. clus- 
ter temperature. We plot data from lLiangetall 1 2000), which were 
partially re-observed by Feretti (2005, in preparation) together with 
data from |Govoni etalj j200l|); |Bacchi etalj <2003|); |Venturi et alj 
Some additional upper limits are collected with the help of 
F. Govoni. We applied a secondary hadronic model as described in 
Dolag & Enfilinl fcOOCh to calculate the radio emission from the sim- 
ulated galaxy clusters. We also added the predictions for the emis- 
sion from primary (upper) and secondary (lower) electrons taken 
from lMiniati et al. (2001). Note that the values for the luminosities 
for primary electrons should be scaled with the electron to proton 
injection ratio R e / P . 
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